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ABSTRACT 


A  mesoscale  investigation  was  conducted  on  the  rapid  coastal  cyclogenesis  that 
occurred  during  Intensive  Observation  Period  (lOP  5  A)  of  the  Experiment  on  Rapidly 
Intensifying  Cyclones  over  the  Atlantic  (ERICA).  Forecasts  from  a  double-nested  version 
of  the  Navy  Operational  Regional  Analysis  and  Prediction  System  (NORAPS)  and 
Multiquadric  Interpolation  (MQI)  objective  analyses  utilizing  operationally  available  and 
some  special  ERICA  data  were  examined  to  study  the  mesoscale  structure  and  frontal 
evolution  associated  with  this  explosively  deepening  coastal  cyclone.  Additionally,  the 
ability  of  NORAPS  to  accurately  simulate  the  explosive  cyclogenesis  was  investigated. 

The  frontal  evolution  showed  characteristics  of  a  classical  occlusion,  similar  to  the 
Norwegian  cyclone  model,  and  marine  frontal  structure  as  described  by  Shapiro  and 
Keyser  (1990).  The  frontal  evolution  was  highly  influenced  by  the  prior  existence  of 
strong  Arctic  and  coastal  fronts.  These  fronts  intensified  during  the  course  of  the  storm 
development  and  did  not  develop  as  a  result  of  the  cyclogenesis. 

The  NORAPS  model  forecasts  were  compared  against  satellite  imagery,  surface 
observations,  MQI  analyses,  and  observed  soundings  taken  during  the  ERICA  study.  The 
double-nested  version  of  NORAPS  was  found  to  be  an  excellent  tool  for  forecasting  the 
mesoscale  frontal  structure  and  intensity  of  this  explosively  deepening  coastal  cyclone. 
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I.  INTRODUCTION 


Rapid  cyclogenesis  at  sea  and  in  coastal  areas  can  pose  a  severe  threat  to  seagoing 
vessels  and  adversely  impact  military  operations.  These  storms  are  associated  with  high 
winds  and  seas  that  also  restrict  visibility  and  pose  substantial  danger  to  both  personnel  and 
property.  Consequently,  accurate  short-range  as  well  as  long-range  forecasts  are  of  great 
importance  to  the  mariner.  With  the  development  of  better  position  and  intensity  forecasts, 
it  is  likely  that  these  storms  will  cause  less  damage  because  better  evasive  and  preparatory 
actions  will  be  possible.  Of  particular  interest  and  concern  are  storms  that  intensify  rapidly. 
Sanders  and  Gyakum  (1980)  define  a  rapidly  deepening  extratropical  cyclone  or  a 
"meteorological  bomb"  as  a  storm  exhibiting  a  central  pressure  fall  exceeding  1  mb  h"’  for 
24  hr  or  more.  While  the  large-scale  processes  that  contribute  to  typical  cyclogenesis  are  well 
imderstood,  it  appears  that  other  processes  beyond  large-scale  baroclinic  dynamics 
contribute  to  the  development  of  rapidly  intensifying  storms.  Although  there  are  variations 
between  individual  storms,  certain  physical  processes  such  as  latent  and  sensible  heat  fluxes 
(Nuss  and  Anthes  1987);  strong  upper-level  forcing  (Uccellini  et  al.  1985;  Wash  et  al.  1988; 
Sanders  1986);  and  a  pre-existing  baroclinic  zone  (Bosart  1981)  have  been  found  to  act 
independently  or  synergistically  to  produce  a  rapid  developer. 

Sanders  (1986)  discussed  the  composite  structure  and  mean  behavior  of  bombs 
forming  in  the  west-central  North  Atlantic  Ocean  during  the  period  January  1 98 1  -  December 
1984.  He  identified  54  cases  of  explosive  deepening  from  the  12  h  operational  analyses 
prepared  by  the  National  Meteorological  Center  (NMC).  Sanders  noted  that  the  average  500 
mb  vorticity  advection  was  highly  correlated  with  the  period  of  maximum  deepening,  which 
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provides  further  evidence  of  the  fundamentally  baroclinic  nature  of  the  bomb.  Sanders 
further  speculated  that  the  large  response  to  the  baroclinic  forcing  is  due  to  the  small  lower- 
tropospheric  static  stability  brought  about  by  surface  latent  and  sensible  heat  fluxes  and  the 
relatively  small  dissipation  over  the  smooth  sea  surface.  Manobianco  (1989)  extended  the 
work  of  Sanders  by  describing  the  three-dimensional  kinematic  and  thermodynamic  synoptic 
and  subsynoptic-scale  structure  for  24  of  the  54  cases  identified  by  Sanders  (1986).  The 
results  of  his  work  support  the  idea  that  explosive  cyclogenesis  is  a  baroclinic  phenomenon 
in  which  the  rapid  development  in  the  presence  of  strong  upper  tropospheric  forcing  appears 
to  be  enhanced  by  a  more  destabilized  lower  troposphere. 

Because  of  inadequate  knowledge  and  parameterizations  of  the  processes  involved, 
operational  models  often  fail  to  adequately  predict  the  development,  intensification  and 
storm  track  of  these  systems  (Bosart  1981).  This  failure  of  numerical  models  to  accurately 
simulate  explosive  cyclogenesis  led  to  the  Office  of  Naval  Research  (ONR)  to  sponsor 
Experiment  on  Rapidly  Intensifying  Cyclones  over  the  Atlantic  (ERICA). 

ERICA  was  conducted  in  the  oceanic  region  east  of  New  England  and  the  middle 
Atlantic  states  and  south  of  Nova  Scotia  and  NeAvfoundland  during  the  period  1  December 
1988  -  28  February  1989.  The  general  area  of  ERICA  measurements  is  shown  in  Fig.  1.  The 
stippled  area  represents  the  region  of  greatest  summed  pressure  falls  for  104  pre-ERICA  type 
storms  and  the  area  that  was  targeted  for  the  highest  concentration  of  ERICA  measurements 
(Hadlock  and  Kreitzberg  1988).  Throughout  this  period,  the  density  of  surface  observations 
over  the  ocean  was  increased  by  measurements  from  specially  deployed  drifting  buoys  and 
ships  of  opportunity.  During  intensive  observation  periods  (lOPs),  additional  low-level  and 
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upper-air  measurements  were  collected  from  instrumented  aircraft,  dropsondes  released 
from  the  aircraft,  and  supplementary  soundings  taken  in  the  eastern  United  States  and 
Canada.  Hadlock  and  Kreitzberg  (1988)  stated  that  the  objectives  of  the  experiment  were 
to: 

(i)  imderstand  the  fundamental  physical  processes  occurring  in  the  atmosphere 
during  rapid  intensification  at  sea; 

(ii)  determine  those  physical  processes  that  need  to  be  incorporated  into 
dynamical  prediction  models  through  efficient  parameterizations  if  necessary 
and; 

(iii)  identify  measurable  precursors  that  might  be  incorporated  in  the  initial 
analysis  for  accurate  and  detailed  operational  model  predictions. 

The  objectives  of  this  paper  are  two  fold:  first,  to  document  the  evolution  and 
mesoscale  features  of  a  rapidly  developing  coastal  cyclone  and  its  associated  frontal 
structures  observed  during  ERICA  lOP  5A  with  all  available  data  and  second,  to  determine 
the  ability  of  the  nested  Navy  Operational  Regional  Atmospheric  Prediction  System 
(NORAPS)  mesoscale  model  to  diagnose  the  fronts  and  cyclone  evolution.  Chapter  II 
describes  the  mesoscale  model  and  analysis  method  used  for  the  study.  Chapter  III  reviews 
the  previous  work  completed  on  lOP  5A  and  describes  model  performance.  Chapter  IV 
contains  the  synoptic  and  mesoscale  description  of  the  event  as  well  as  model  verification 
through  actual  observations.  A  summary  and  conclusions  is  found  in  Chapter  V. 
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11.  MESOSCALE  MODEL  AND  ANALYSES 


A.  MESOSCALE  FORECAST  MODEL 

NORAPS  version  6.1  is  a  globally  relocatable,  triple-nested  mesoscale  model  run 
operationally  by  the  Naval  Fleet  Numerical  Meteorology  and  Oceanography  Center  for 
several  areas  of  the  world;  the  continental  U.S.  (including  the  Western  Atlantic  Ocean); 
Indian  Ocean;  Europe  (including  the  Mediterranean);  and  Asia  ( including  the  Western 
Pacific  Ocean).  The  NORAPS  model  has  been  described  by  Hodur  (1987)  and  more  recently 
by  Liou  et  al.  (1994).  Initial  conditions  (analyses)  are  prepared  on  16  standard  pressure 
levels  through  multivariate  optimum  interpolation  (MVOI)  (Barker  1992).  Model  forecasts 
are  obtained  by  integrating  the  hydrostatic  primitive  equations  and  physical  parameterizations 
on  a  staggered  Arakawa  C-grid.  The  model  atmosphere  is  divided  into  24  layers  from  the 
surface  to  50  mb.  Of  the  24  levels,  7  are  concentrated  below  850  mb  to  provide  adequate 
resolution  of  the  planetary  boundary  layer  (PBL).  Physical  parameterizations  for  NORAPS 
include:  radiative  transfer,  Kuo  (1974)  type  cumulus  parameterization  for  deep  convection, 
and  Tiedtke  et  al.  (1989)  parameterization  for  large-scale  condensation.  The  multi-level 
boundary  layer  is  parameterized  by  using  K-theory  for  the  surface  (Louis  1979)  and  a  1.5 
turbulent  kinetic  energy  closure  scheme  for  vertical  eddy  flux  parameterization  (Langland 
and  Liou  1994).  Precipitation  falling  into  unsaturated  layers  is  partially  evaporated, 
depending  on  the  relative  humidity  of  the  subcloud  layers.  NORAPS  grids  can  be  selected 
using  any  of  the  conformal  projections:  Lambert,  Mercator,  or  polar  stereographic.  At  inner 
lateral  boundary  zones,  terrain  heights  of  the  inner  fine  mesh  are  matched  with  the  terrain 
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heights  of  the  coarser  mesh  so  that  sigma  coordinates  are  consistent  in  those  regions.  During 
time  integration,  the  boundary  conditions  for  the  two  inner  meshes  are  updated  with  each 
timestep,  while  the  coarser  outer  grid  is  updated  every  6  or  12  h  by  the  forecasts  of  the  Navy 
Global  Atmospheric  Predictions  System  (NOGAPS)  (Liou  et  al.  1994).  The  U.  S.  Navy's  ten 
minute  data  base  provides  the  terrain  fields  used  by  NORAPS 

A  double-mesh  version  of  NORAPS  has  been  used  in  this  study.  The  forecast 
domain  for  the  fine  mesh  contains  121  x  121  grid  points  with  a  grid  size  of  20  km  centered 
at  45  °N,  65  °  W.  The  domain  for  the  coarse  mesh  contains  109  x  89  grid  points  with  a  grid 
size  of  60  km  centered  at  45  °  N,  65  °  W  (Fig.  2).  The  forecast  model  was  initialized  with  60 
km  NORAPS  01  analyses.  The  initialization  was  supplemented  by  regional  operational 
surface  data,  rawinsonde  data,  and  special  ERICA  data.  The  model  simulation  started  at  1200 
UTC  20  January  and  was  integrated  36  h  to  0000  UTC  22  January  1989. 

B.  SURFACE  ANALYSES 

In  order  to  verify  the  model  forecasts,  reanalyses  of  the  sea  level  pressure  and  surface 
temperatures  were  produced  using  a  multiquadric  interpolation  (MQI)  scheme  on  the 
NORAPS  60  km  and  20  km  grids.  As  described  by  Nuss  and  Titley  (1994),  the  MQI  fits 
surfaces  through  all  observations  in  such  a  way  that  values  can  be  interpolated  to  any  point 
on  a  grid.  NORAPS  analysis  fields  were  used  as  a  first-guess.  The  degree  of  surface  fit  to  the 
observations  is  determined  by  specification  of  filtering  and  smoothing  parameters.  The  user 
may  vary  the  root  mean  square  difference  between  Observations  and  the  first-guess  field. 
Similarly,  the  smoothing  of  the  observations  may  be  varied  in  the  analysis.  Nuss  and  Titley 
have  demonstrated  that  the  MQI  technique  has  skill  in  the  interpolation  of  scattered  data. 
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in.  REVIEW  AND  METHODOLOGY 


The  time  period  18-22  January  1989  was  one  of  the  most  active  during  the  ERICA 
project.  During  this  five  day  period,  two  rapidly  intensifying  cyclones  developed  within  36 
h  over  the  western  North  Atlantic  Ocean.  The  first  cyclone,  lOP  5,  developed  off  the  mid- 
Atlantic  coast  and  tracked  along  the  north  wall  of  the  Gulf  Stream,  reaching  a  maximum 
deepening  rate  of  36  mb  (18  h)’*  (Hadlock  et  al.  1989).  The  second  cyclone,  lOP  5A, 
developed  northeast  of  Lake  Ontario  and  tracked  along  the  Canadian  Maritimes  coast, 
reaching  a  maximum  deepening  rate  of  21  mb  (9  h)‘*  (Fig.  3).  This  latter  storm  is  of  particular 
interest  because  it  is  a  coastal  developing  storm,  which  exhibited  characteristics  of  rapid 
deepening  over  land  as  described  by  Mass  and  Schultz  (1993),  and  marine  development  as 
summarized  by  Shapiro  and  Keyser  (1990).  Of  particular  interest,  pre-existing  frontal  zones 
intensified  during  the  rapid  cyclogenesis  and  evolved  into  cyclone  fi-ontal  structures  rather 
than  forming  as  a  result  of  the  cyclogenesis. 

A.  REVIEW  OF  PREVIOUS  WORK 

Spinelli  (1992)  conducted  a  synoptic  investigation  of  lOP  5  A  and  determined  several 
factors  that  contributed  to  the  rapid  development  of  the  cyclone.  Those  factors  were:  (i) 
significant  lower  tropospheric  thermal  advection  preceding  the  rapid  intensification;  (ii) 
favorable  superposition  of  a  mobile  500  mb  trough  over  the  fi’ontal  wave  providing  upper- 
level  support;  (iii)  the  presence  of  a  jet  streak  on  the  eastern  side  of  the  300  mb  trough  and; 
(iv)  intense  upward  vertical  motions  within  the  frontal  cloud  band  of  the  cyclone.  In  addition, 
an  evaluation  of  the  60  km  NORAPS  model  forecast  versus  the  60  km  NORAPS  01 
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Figure  3:  NMC  SLP  (mb)  vs  time  for  (a)  ERICA  lOP  5  (b)  ERICA  lOP  5A 
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analyses  was  conducted  to  determine  the  performance  of  the  model  in  simulating  a  rapidly 
developing  coastal  storm. 

In  evaluating  the  NORAPS  model,  Spinelli  (1992)  noted  that  the  NORAPS  60  km 
models  forecasts  of  the  ^optic-scale  features,  (e.g.,  jet  streaks,  vorticity  centers,  upper- 
level  troughs)  position  and  intensity  were  “relatively  accurate.”  However,  the  60  km 
NORAPS  analyses  and  forecasts  were  too  coarse  to  resolve  the  mesoscale  features  that  may 
have  played  an  important  role  in  the  rapid  development  of  the  lOP  5  A  cyclone.  To  resolve 
these  features,  Spinelli  performed  a  hourly  subjective  analysis,  which  revealed  the 
development  of  a  secondary  low  pressure  center  well  to  the  southeast  of  the  primary  low  (Fig. 
4).  This  secondary  development  was  not  indicated  in  the  synoptic  NORAPS  andyses.  Based 
on  the  topography  of  the  region  (Fig.  5),  Spinelli  speculated  that  this-low  was  probably  a 
result  of  lee  cyclogenesis  on  the  lee  (eastern)  side  of  the  northern  Appalachian  Mountains. 
The  coastal  low  continued  to  progress  northeast  and  deepen  rapidly  and  likely  became  the 
primary  low  center  for  the  lOP  5A  cyclone.  Figure  6  shows  the  tracks  of  the  two  separate 
lows  from  20/1 500Z  to  21/0300Z  as  derived  from  the  subjective  analyses.  The  6-h  positions 
of  the  lOP  5  A  cyclone  as  analyzed  by  the  60  km  NORAPS  are  represented  by  the  open  square 
symbols.  One  significant  result  of  these  detailed  analyses  is  that  there  is  a  clear  indication  that 
mesoscale  processes  not  necessarily  captured  by  the  60  km  NORAPS  analysis  or  forecasts 
were  occurring  during  the  time  period  examined. 

Further  assessment  of  the  NORAPS  model  to  predict  mesoscale  features  was  also 
conducted  by  Cameron  (1993).  Specifically,  a  comparison  between  the  frontal  evolutions 
associated  with  other  rapidly  deepening  storms  and  the  coastal  lOP  5  A  storm  was  performed 
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Figure  4.  Subjective  surface  pressure  analysis  (solid,  contour  interval  4  mb)  with 
subjective  temperature  analysis  added  (dashed,  contour  interval  2°C)  at  20/1 5Z  January 
1989  (from  Spinelli  1992). 


to  determine  if  maritime  frontal  evolution  can  be  generalized  to  a  coastal  environment.  In 
his  investigation,  Cameron  found  that  the  frontal  evolution  of  lOP  5  A  did  not  neatly  fit  into 
a  angle  category  of  development.  Rather,  the  storm  exhibited  properties  similar  to  a  classical 
occlusion  and  some  evidence  of  the  bent-back  warm  front  structure  described  by  Shapiro  and 
Keyser  (1990).  For  example,  composite  of  the  Arctic  front  evolution  as  depicted  by  the  60 
km  NORAPS  model  is  shown  in  Fig.  7  at  12  h  intervals  from  20/1200Z  to  22/OOOOZ.  Figure 
8  is  a  coastal  front  composite  every  ax  hours  from  20/1200Z  to  22/OOOOZ.  These  composites 
show  that  strong  Arctic  and  coastal  fronts  present  at  the  incipient  stage  of  cyclone 
development  intensified  as  a  result  of  the  rapid  cyclogenesis  and  became  an  integral  part  of 
the  synoptic-scale  frontal  structure. 

B.  NESTED  NORAPS  MODEL  PERFORMANCE 

A  significant  result  of  these  two  previous  studies  was  that  the  need  for  a  higher 
resolution  mesoscale  model  to  accurately  predict  secondary  cyclogenesis  and  fiirther 
investigate  the  mesoscale  features.  This  study  utilizes  the  results  from  a  36-h  NORAPS 
simulation  on  the  20  km  mesh  previously  described.  A  20  km  MQI  analysis  of  the  SLP  and 
surfece  temperatures  is  also  utilized  to:  determine  the  SLP  and  positions  of  lOP  5  A;  perform 
frontal  analysis;  and  to  validate  NORAPS  output.  In  evaluating  forecast  model  performance, 
time  series  from  Sable  Island,  MQI  and  National  Meteorological  Center  (NMC)  analyses,  in 
situ  observations,  and  the  NORAPS  60  km  analysis  were  used.  In  this  manner,  the  model  skill 
in  predicting  the  evolution  and  mesoscale  features  of  the  lOP  5  A  cyclone  wiU  be  compared 
closely  to  the  observations. 
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Figure  6.  Tracks  (solid  lines  connecting  open  circles)  and  central  pressures  (mb)  of 
primary  and  secondary  lows  from  20/1500Z  to  21/0300Z.  The  hour  and  last  two  digits  of 
the  central  pressure  are  shown  above  and  below  the  position  markers.  Open  squares  are 
positions  of  main  low  derived  from  NORAPS  analyses. 


1.  Position  and  Intensity 

NORAPS  did  an  excellent  job  in  forecasting  the  rapid  development  of  the  lOP  5  A 
cyclone.  A  comparison  of  intensity  and  position  forecasts  is  provided  in  Table  1.  While  the  6-h 
forecast  central  pressure  of  the  low  was  2  mb  too  high,  the  model  was  2  mb  too  deep  by 
tau  24,  but  remained  2  mb  weak  for  the  remaining  period.  The  36-h  analysis  and  NORAPS 
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Figure  7.  Composite  Arctic  front  positions  in  twelve  hour  increments  from  20/12Z  to 
22/OOZ  January  1989  (from  Cameron  1993). 


90  80  70  60  50 

Figure  8.  Composite  coastal  front  positions  in  six  hour  increments  from  20/1 2Z  to 
22/OOZ  January  1989  (from  Cameron  1993). 
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forecast  positions  of  the  ERICA  lOP  5  A  cyclone  at  6-h  intervals  starting  from  20/1 200Z  are 
depicted  in  Fig.  9.  The  initial  analyzed  position  is  just  north  of  Lake  Ontario.  Note  that  the 
Table  1.  Position  and  intensity  forecasts  for  lOP  5 A.  Position  errors  are  read  as  the  distance 
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minimum  central  pressure  of  964  mb  occurs  at  21/1 800Z  close  to  the  southwestern  edge  of 
Newfoundland.  The  storm  deepens  31  mb  in  the  18-h  period  from  21/OOOOZ  to  21/1800Z 
with  a  model  pressure  difference  of  +2  mb  (Fig.  10).  After  18  h  (21/0600Z),  the  forecast 
shows  an  increasing  bias  in  the  Qrclone  position,  with  NORAPS  tending  to  place  the  storm 
center  somewhat  to  the  northwest  of  the  observed  position. 

2.  Time  Series  Comparisons 

To  further  examine  model  performance,  a  time  series  of  SLP,  surface  temperature, 
wind  speed  and  direction  from  Sable  Island  (WSA)  during  a  cold  frontal  passage  will  be 
examined.  This  station  is  located  on  a  small  island  (43.9°  N,  60.0°W),  while  the  available 
every  hour,  and  are  plotted  with  the  hourly  model  output  in  Fig.  1 1 .  A  point  verification  is 
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a  challenging  test  for  a  numerical  weather  prediction  model. 

In  comparing  the  time  series  in  Figs.  1 1  a-b,  NORAPS  generally  has  similar  SLP  but 
warmer  near-surface  temperatures  when  compared  to  the  observations  until  frontal  passage, 
which  occurs  22  h  (21/1 OOOZ)  into  the  simulation.  The  surface  observations  show  a  sharp 
jump  in  temperature  near  the  frontal  passage,  which  is  captured  skillfully  by  the  model.  Prior 
to  frontal  passage  (21/OOOOZ  -  21/0900Z),  NORAPS  shows  a  24  mb  (9  h)‘*  pressure  fall  as 
compared  to  the  observational  rate  of  27  mb  (9  h)'^  After  frontal  passage,  the  model 
remains  approximately  5  mb  too  high  at  this  point  for  the  remainder  of  the  period.  The  earlier 
warm  bias  of  the  suifece  temperatures  begins  to  diminish  and  NORAPS  forecasts  the  surface 
temperatures  after  frontal  passage  exceedingly  well. 

NORAPS  captures  the  observed  wind  speed  and  direction  well  (Figs.  1 1  c-d).  The 
NORAPS  wind  speed  is  nearly  equal  with  the  observations  until  frontal  passage.  At  frontal 
passage,  the  observations  show  a  momentary  drop  in  wind  speed  with  a  gradual  veering  of 
the  winds  to  the  west.  NORAPS  follows  both  of  these  trends  except  that  the  wind  speed 
remains  nearly  steady  during  frontal  passage  and  remains  somewhat  weak  for  the  remaining 
period. 

C.  SUMMARY  OF  MODEL  PERFORMANCE 

From  previous  investigations  of  the  lOP  5A  cyclone,  it  was  demonstrated  that  the 
coarse  grid  NORAPS  model  produced  a  realistic  simulation  of  the  rapid  cyclogenesis  event. 
The  important  large-scale  features  such  as  the  location  and  intensity  of  the  500  mb  and  300 
mb  troughs,  vorticity  centers,  jet  streaks  and  upper-level  divergence  areas  were  all 
represented  well.  The  reader  is  referred  to  Spinelli  (1992)  and  Cameron  (1993)  for  a  detailed 
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Figure  9.  Sea  surface  temperatures  (solid,  TC)  on  20  January  1989.  NORAPS  (crosses) 
and  the  MQI  analysis  (circles)  6-h  surface  cyclone  positions  and  central  sea  level  pressures 
between  12  UTC  20  January  and  0000  UTC  22  January  1989. 


SLP  VS  Time 


Figure  10.  Analyzed  (MQI)  and  simulated  sea  level  pressure  (mb)  at  the  center  of  the 
cyclone  for  the  36-h  simulation  period. 
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Temperature  VS  Time  (WSA) 


SLP  VS  Time  (WSA) 


Wind  Speed  VS  Time  (WSA) 


Wind  Direction  VS  Time  (WSA) 


Figure  11.  Time  series  comparisons  from  Sable  Island  (WSA)  during  cold  frontal 
passage:  (a)  sea  level  pressure  (mb);  (b),  temperature  (TC);  (c),  wind  speed  (m/s);  (d), 
wind  direction  (degrees). 
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analysis  of  the  coarse  grid  NORAPS  model  performance  in  simulating  ERICA  lOP  5  A. 

Although  it  is  evident  that  the  overall  performance  of  the  coarse  grid  NORAPS  model 
in  predicting  the  rapid  cyclogenesis  was  successhil,  a  higher  resolution  model  was  used  to 
investigate  the  mesoscale  features  of  lOP  5A.  In  comparing  the  forecast  fields  against 
observations,  it  is  shown  that  the  latest  version  of  NORAPS  demonstrated  considerable  skill 
in  forecasting  the  mesoscale  features  of  lOP  5 A.  Moreover,  the  model  simulation  provides 
a  high  resolution  and  dynamically  consistent  data  set.  Therefore,  it  will  be  used  in  the 
subsequent  analyses  of  the  structure,  development  and  evolution  of  lOP  5  A. 


22 


IV.  MESOSCALE  EVOLUTION 


The  following  sections  present  the  lOP  5A  cyclone  development  and  frontal 
evolution  in  three  stages:  the  incipient  development  stage  (20/1200Z  -  21/OOOOZ);  the 
explosive  development  stage  (21/OOOOZ  -  21/1800Z)  and  the  mature  cyclone  stage 
(21/1800Z-22/0000Z),  Six  hour  surface  temperature  and  pressure  analyses  (MQI)  are 
presented  in  Fig.  12  along  with  the  NORAPS  6  h  forecasts  for  the  same  periods.  Frontal 
evolution  will  be  described  to  address  specific  details  of  lOP  5A’s  fi-ontal  life  cycle  from  its 
incipient  phase  through  warm-core  seclusion.  The  frontal  analysis  is  based  on  a  detailed 
examination  of  the  surface  observations  obtained  from  the  ERICA  dataset. 

The  rate  of  change  of  the  gradient  of  potential  temperature  following  a  parcel  is  used 
to  examine  the  kinematics  of  the  simulated  frontogenesis.  The  horizontal  frontogenesis 
function  in  pressure  coordinates  is  defined  as: 


£ 

dt 


^  ^  dx  dx 


^.d6y5v.._.-d6  89  .dv 
dy  dy  dx  dy  dx  dy 


dx  dp  dx  dy  dp  dy  dx  dx  dt  dy  dy  dt 


(1) 


For  this  study  only,  the  confluence  and  shear  terms,  the  first  two  terms  on  the  right  hanH  side 
of  the  equation,  will  be  used  to  estimate  d/dt|  Vp0 1 .  The  remaining  terms  were  not  found  to 
be  significant.  The  1000  mb  level  of  the  NORAPS  model  was  used  to  produce  fields  of 
surface  frontogenesis.  In  addition,  model  vertical  soundings  will  be  evaluated  and  compared 
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to  actual  soundings  observed  during  the  experiment.  Model  cross  sections  will  also  aid  in 
determining  frontal  position  and  structure. 

There  remains  significant  interest  in  frontal  evolution  as  noted  in  Neiman  and  Shapiro 
(1993),  Mass  and  Schultz  (1993),  and  McGinnigle  (1988).  In  addition,  forecasting  the 
structure  and  development  of  sub-synoptic  scale  features  in  coastal  areas  remains  a  difficult 
task,  largely  because  timely  synoptic  observations  end  at  the  shoreline.  Even  with  special 
collection  efforts  such  as  the  ERICA  study,  large  ocean  areas  or  sparsely  populated  land  areas 
may  not  be  adequately  sampled,  leaving  large  gaps  where  significant  mesoscale  activity  may 
develop.  In  an  effort  to  combat  the  sparsity  of  data,  gridded  data  sets  are  often  used  in  place 
of  unverified  analyses.  Therefore,  this  study  utilizes  a  data  base  consisting  of  hourly 
NORAPS  forecast  fields  augmented  with  the  special  data  available  from  the  ERICA  studies. 
The  ability  of  the  high  resolution  model  to  resolve  the  mesoscale  features  and  frontal 
evolution  of  a  rapidly  deepening  cyclone  will  continue  to  be  documented  with  detailed  model 
station  soundings  and  cross  sections  compared  with  the  observed  fields. 

A.  INCIPIENT  DEVELOPMENT  STAGE  (2O/12OOZ-21/0000Z) 

Horizontal  surface  analysis  of  pressure  and  temperature  with  associated  frontal 
positions  are  presented  in  6-h  intervals  to  illustrate  the  cyclone's  mesoscale  frontal  stmcture 
during  this  period  (Figs.  12a,c,e,g,i,k,m)  (left  panels).  NORAPS  simulated  fronts  are  also 
depicted  (Figs.  12b,d,f,hj,l,n)  (right  panels)  for  comparison.  The  36-h  analyzed  and  forecast 
storm  track  shown  in  Fig.  9  depicts  sea  surface  temperatiores  valid  at  initialization.  Of 
particular  note  is  the  strong  temperature  gradient  marking  the  Gulf  Stream  off  the  coast  of 
North  Carolina.  The  ice  edge  is  also  present  further  north  along  the  coast  of  eastern  Quebec 
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Figure  12.  Six-hourly  surface  MQI  analyses  (left  panels)  and  NORAPS  forecasts  (right 
panels)  of  pressure  (mb,  contour  interval  2  mb)  and  temperature  (dashed,  contour  interval 
2°  C)  for  20  January  1989;  (a-b)  and  (c-d),  1200Z  and  1800Z  respectively,  21  January; 

(g"h)>  (i-j)?  and  (k-1)  OOOOZ,  0600Z,  1 200Z,  and  1 800Z  respectively,  (m-n) 
22/OOOOZ  January  1989.  Fronts  (conventional  symbols)  are  also  depicted. 
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22/OOZ  sft  sip  MQI _  22/OOZ  sft  sip  NORAPS 


and  continuing  south  to  the  tip  of  Newfoxmdland.  Not  depicted  is  the  broken  ice  coverage  in 
the  northern  portion  of  the  Gulf  of  St.  Lawrence  near  the  coast  of  southern  Quebec. 

At  the  beginning  of  the  period  (Fig.  12a),  the  surface  low  is  located  just  north  of  Lake 
Ontario  with  a  central  pressure  of  1 003  mb.  Further  to  the  east,  the  lOP  5  storm  is  weakening 
southeast  of  Newfoundland.  A  complex  pattern  with  four  distinct  baroclinic  zones  are 
evident.  In  addition  to  the  primary  warm  and  cold  fronts  associated  with  the  incipient  lOP 
5A  cyclone,  a  pre-existing  Arctic  front  extends  southwest  from  central  Quebec,  upstream 
of  the  low  through  Michigan,  Illinois  and  ending  in  northeastern  Missouri  (full  extent  not 
shown).  This  front  is  very  intense  north  of  the  Great  Lakes  as  evidenced  by  the  tight 
temperature  gradient.  A  weak  northern  warm  frontal  band  extends  from  the  low  center,  east 
through  Quebec  and  into  New  Brunswick.  The  polar  front  extends  south  of  the  low  along  the 
pressure  trough  through  western  New  York,  central  Pennsylvania,  Virginia  and  across 
western  North  Carolina.  This  cold  front  is  not  clearly  seen  in  the  surface  temperature  analysis 
at  this  time,  but  is  evident  in  the  pressure  pattern  and  is  carried  on  operational  NMC  charts 
(not  shown).  Additionally,  a  pre-existing  mesoscale  coastal  front  off  the  North  Carolina  coast 
associated  with  strong  sea  surface  temperature  gradients  is  evident.  A  very  weak  warm  front 
intersects  the  coastal  front  along  the  temperature  ridge  line  and  extends  to  the  southeast. 
There  is  also  a  weak  coastal  front  beginning  to  form  along  the  southern  coast  of  Maine, 
which  extends  to  Massachusetts.  There  is  a  warm  tongue  along  the  Southeast  Coast.  This 
feature,  along  with  the  Arctic  and  polar  fronts  are  evident  at  850  mb  as  well  (Fig.  13). 

By  20/1 800Z,  the  low  center  moves  east  and  deepens  five  mb  to  998  mb  (Fig.  12c). 
The  northern  portion  of  the  Arctic  front  moves  east  at  20  kts  and  catches  up  with  the  surface 
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Figure  13.  850  mb  temperature  (solid,  positive  values,  dashed,  negative  values,  contour 
interval  1  °C)  valid  at  20/1200Z  January  1989. 

low  in  southern  Quebec.  The  portion  of  the  front  southwest  of  Lake  Ontario  began  to 
dissipate.  A  weak  warm  front  is  still  evident  north  and  east  of  the  surface  low.  This  front 
intersects  the  Arctic  front  at  the  low  position  in  southern  Quebec  and  extends  east,  crossing 
New  Brunswick  into  the  partially  ice  covered  Gulf  of  St.  Lawrence.  The  polar  front  pushes 
east  and  extends  southwest  into  South  Carolina.  The  thermal  packing  of  this  front  is  more 
noticeable  at  this  time.  The  coastal  front  and  associated  warm  front  are  notably  stronger  as 
well.  Figure  14  depicts  surface  frontogenesis  at  6  h  intervals  starting  at  20/1 200Z  and  ending 
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Figure  14.  (Continued) 
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at  22/OOOOZ.  Figure  14b  shows  five  areas  of  positive  values.  These  areas  correlate  directly 
with  the  locations  of  the  Arctic  fi'ont;  polar  front,  the  associated  warm  front;  the  two  coastal 
fronts  off  the  Virginia  coast  and  southern  Maine;  and  the  strongly  baroclinic  ice  edge  off 
eastern  Canada  to  the  north.  The  NORAPS  depicts  the  location  of  each  of  these  fronts  well. 
Of  particular  note  is  how  well  NORAPS  represents  the  position  and  thermal  packing  behind 
the  polar  and  coastal  fronts.  The  analysis  puts  the  polar  front  in  the  pressure  trough  just 
ahead  of  the  thermal  packing  located  in  southeastern  Virginia.  Further  east  there  is  a 
weakening  of  the  thermal  gradient  just  off  the  coast  and  then  the  thermal  gradient  strengthens 
behind  the  coastal  front.  NORAPS  shows  these  features  well.  Satellite  imagery  (Fig.  15)  at 
20/1 80 IZ  is  in  agreement  with  the  position  of  both  fronts  and  shows  the  distinct  separation 
of  the  comma  cloud  mass  associated  with  the  incipient  low  and  polar  front  to  the  north  and 
the  cloud  mass  to  the  southeast  associated  with  the  coastal  front. 

By  21/OOOOZ,  the  synoptic-scale  low  moves  to  the  southeast  and  is  positioned  over 
central  Maine  with  a  pressure  of  995  mb.  The  mesoscale  details  associated  with  this 
movement  will  be  discussed  shortly.  Rapid  deepening  of  the  storm  is  now  beginning  to 
occur.  The  Arctic  front  has  moved  east  as  well  (Fig.  12e)  and  is  located  just  west  of  the  low. 
Over  the  last  6  h,  the  polar  front  has  weakened  considerably  and  has  combined  with  the 
coastal  front,  which  has  intensified  over  the  western  North  Atlantic  to  form  an  occlusion. 
This  occlusion  is  classified  as  a  warm  occlusion  since  the  coldest  air  is  ahead  of  the  warm 
front.  This  newly  combined  frontal  system  will  be  referred  to  as  the  synoptic  cold  front 
hereafter.  The  synoptic  warm  frontal  boundary  associated  with  this  cold  front  can  clearly  be 
seen  in  the  isotherm  field  and  on  the  surface  frontogenesis  chart  as  a  separate  area  of 
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Figure  15.  Goes  enhanced  IR  satellite  imagery  at  20/1 80 IZ  January  1989. 

positive  values  (Fig.  14c).  The  surface  frontogenesis  also  indicates  that  the  previous  synoptic 
warm  frontal  band  associated  with  the  incipient  low  has  now  moved  along  the  intensifying 
baroclinic  zone  along  the  northern  Gulf  of  St.  Lawrence  and  southern  Quebec  and  become 
stationary.  This  increased  intensification  of  the  baroclinic  zone  in  the  northern  portion  of  the 
Gulf  of  St.  Lawrence  is  due  to  greater  convergence  along  the  coastal  mountain  ridge 
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associated  with  the  eastward  moving  cyclone. 

An  examination  of  the  2 1/OOOOZ  850  mb  temperatures  (Fig.  1 6)  reveals  that  a  distinct 
rotation  of  the  thermal  ridge  over  New  England  from  a  southwest-northeast  orientation  to  a 
more  north-south  alignment  had  occurred  over  the  last  6  h.  The  presence  of  the  warm 
tongue  has  become  more  pronounced  as  can  be  seen  by  the  6°  C  contour  line  which  now 
extends  north  of  40  °N. 

GOES  infrared  satellite  imagery  valid  at  21/OOOlZ  (Fig.  17)  shows  a  distinct  comma 
head  with  the  largest  cloud  mass  ahead  of  the  low,  over  New  Brunswick  and  southeast  of 
Novia  Scotia.  A  developing  dry  slot  south  of  the  cyclone  is  also  evident.  This  mature  comma 
cloud  structure  formed  vviien  the  comma  cloud  pattern  in  the  polar  air  upstream  of  the  main 
frontal  band  merged  with  the  main  cloud  mass  associated  with  the  old  coastal  front.  The 
cloud  mass  merging  and  evolution  during  the  early  part  of  the  rapid  deepening  is  similar  to 
the  "instant  occlusion"  scenario  as  described  by  Mullen  (1983).  The  cloud  band  in  Fig.  17 
shows  a  distinct  western  border  marking  the  limiting  streamline  of  the  warm  conveyor  belt 
structure  as  described  by  Carlson  (1980),  which  separates  the  dry  and  moist  air  streams. 

The  NORAPS  forecasts  of  the  temperature  and  wind  in  the  vertical  will  now  be 
assessed.  To  do  this,  a  vertical  cross  section  (Fig.  18)  at  21/OOOZ  from  Brookhaven,  Long 
Island  (BNL)  40.9°  N,  72.9°W  to  Sable  Island  (WSA)  43.9°  N,  60.0°W  is  compared  to  the 
model  cross  section  in  Fig.  19.  This  section  contains  five  stations  and  provides  adequate 
resolution  in  the  horizontal.  The  location  of  the  cross  section  is  indicated  in  Fig.  12e. 

The  left  edge  of  the  cross  section  in  Fig.  18  depicts  the  weak  coastal  baroclinic  zone 
southwest  of  the  coastal  low.  The  cross  section  depicts  the  occluded  front  near  station  N234. 
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12  21/OOZ  NORAPS  850  mb  temps 


Figure  16.  850  mb  temperature  as  in  Fig.  13,  except  for  21/OOOOZ  January  1989. 


The  nearly  vertical  theta  lines  in  the  western  part  of  the  cross  section  represent  the  leading 
edge  of  continental  polar  air  that  has  now  moved  southward  from  eastern  Canada.  The 
location  of  the  occlusion  is  clearly  evident  in  the  tongue  of  warm  air  near  station  N234 
approximately  15  nm  south  of  Maine.  Further  east,  the  cross  section  terminates  in  the  lower 
theta  air  just  north  of  the  warm  front.  NORAPS  moisture  values  indicate  80%-90%  saturation 
in  the  vicinity  of  the  warm  tongue  while  the  observed  soundings  only  indicate  50%-80%. 
The  model  cross  section  (Fig.  19)  describes  all  of  these  features  generally  well  and  of 
particular  note  is  how  well  NORAPS  handles  the  change  in  wind  direction  to  the  south  at  the 
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Figure  17.  Goes  enhanced  IR  satellite  imagery  at  21/OOOlZ  January  1989. 
base  of  the  warm  tongue. 

The  NORAPS  forecast  fields  and  satellite  data  illustrate  the  complexities  and 
uniqueness  of  the  lOP  5A  cyclone.  Specifically,  this  cyclone  does  not  form  along  a  single, 
broad  baroclinic  zone  as  discussed  by  Neiman  and  Shapiro  (1993)  or  by  Mass  and  Schultz 
(1993).  In  particular,  the  cloud  and  frontal  analysis  show  several  important  pre-existing 
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baroclinic  zones  near  llie  incipient  cyclone,  an  intensifying  coastal  front  with  a  distinct  warm 
tongue  to  the  south  and  a  strong  Arctic  front  present  west  of  the  low. 

B.  EARLY  MESOSCALE  DEVELOPMENT 

Spinelli  (1992)  performed  a  subjective  hourly  hand  analysis  which  revealed  the 
presence  of  a  secondaiy  low  pressure  center  well  to  the  southeast  of  the  primary  center  just 
prior  to  explosive  deepening  (20/1 500Z)  (Fig  4).  Her  analysis  of  this  secondary  low  provided 
the  motivation  for  further  analysis  of  this  area  using  an  objective  technique.  For  this  study, 
an  objective  analysis  at  3  h  intervals  was  performed  on  the  data  using  the  Multiquadric 
Interpolation  technique  as  described  by  Nuss  and  Titley  (1994).  A  mesoscale  comparison  of 
how  well  NORAPS  forecast  fields  compared  to  MQI  analyses  was  performed.  Comparisons 
of  the  MQI  analysis  and  the  NORAPS  model  output  with  verifying  observations  are  depicted 
in  Fig.  20  for  the  period  of  20/1200Z  to  21/0300Z.  The  following  discussion  will  focus  on 
the  MQI  analysis  (left  panels)  with  comments  on  the  NORAPS  model  (right  panels)  where 
they  differ  significantly. 

From  20/1200Z  to  20/1 800Z  (Figs.  20a,  c,e)  the  primary  lOP  5 A  low  quickly  tracked 
to  the  northeast  at  a  speed  of  35  kts  and  deepened  from  1003  mb  to  998  mb.  At  the  same 
time,  a  developing  low  pressure  area  in  southeastern  New  York  also  moved  northeast  and 
was  located  near  42.7N  72.5W  (southern  Vermont)  at  20/1 800Z.  While  the  wind  pattern  did 
not  indicate  a  cyclonic  circulation,  minimum  pressure  in  the  area  dropped  from 
approximately  1007  mb  to  1002  mb  during  the  3  h  period.  The  NORAPS  depiction  agrees 
well  with  the  analysis  through  20/1 800Z. 

From  20/1 800Z  to  20/2100Z  the  cyclone  moved  southeast  to  45.5  °N,  71 .4°  and 
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valid  at  21/1200Z  January  1989. 
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remained  at  998  mb.  The  low  pressure  area  to  the  south  moved  eastward  to  near  the 
southwest  comer  of  Maine  and  SLP  dropped  in  the  area  down  to  999  mb  (Fig.  20g).  This 
represents  a  deepening  of  9  mb  6  h’’  for  this  coastal  region.  The  observations  now  begin  to 
indicate  a  cyclonic  circulation,  showing  a  distinct  wind  shift  fi'om  northerly  at  Manchester 
Airpark,  NH  (MHT)  to  westerly  at  Beverly  Municipal,  MA  (B  VY)  to  southerly  at  the  buoy 
lOSN  with  lowest  pressures  being  reported  by  the  coastal  stations  (station  identification  not 
shown). 

The  fine  mesh  NORAPS  model  did  an  exceptional  job  at  forecasting  the  secondary 
cyclogenesis.  The  corresponding  NORAPS  depiction  of  the  secondary  development  is  best 
seen  in  Figs.  20  f,h.  The  area  of  broad  troughing  is  very  similar  to  the  analyzed  troughing  6 
h  into  the  simulation.  By  20/2  l  OOOZ,  NORAPS  develops  a  1001  mb  secondary  cyclone  near 
43. 3N  70.7W  (southwest  Maine),  while  the  MQI  analysis  is  showing  a  strong  pressure  trough 
with  approximately  999  mb  for  the  area.  Although  2  mb  weak  in  pressure,  NORAPS  does 
an  excellent  job  in  positioning  the  new  low  just  25  nm  northeast  of  the  cyclonic  circulation 
evident  in  the  observations  previously  described.  Although  this  low  is  not  specifically 
indicated  in  the  analysis,  the  analyzed  pressures  for  this  area  are  more  representative  of  the 
observations.  Due  to  the  univariate  nature  of  the  MQI  analysis,  it  may  of  not  been  capable 
to  place  a  low  center  there  based  only  on  the  pressure  field.  NORAPS  continues  to  move  the 
northern  primaiy  cyclone  to  the  northeast  and  deepens  it  to  998  mb.  Based  on  the  topography 
(Fig.21),  and  the  sub-synoptic  scale  flow  in  the  area,  it  is  likely  that  this  secondary  low 
pressure  area  is  a  result  of  lee  cyclogenesis  on  the  lee  side  of  the  northern  Appalachian 
Moimtains. 
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Figure  20.  Three-hourly  surface  MQI  analyses  (left  panels)  and  NORAPS  forecasts  (right 
panels)  (isobars,  solid,  contour  interval  1  mb)  with  verifying  observations  20  Januaiy  1989;  (a- 
b),  (c-d),  (e-f),  and  (g-h)  1200Z,  1500Z,  1800Z,  and  2100Z  respectively,  21  January;  (i-j),  (k-1) 
(m-n),  OOOOZ,  0300Z,  and  0600Z,  respectively.  Six-hour  intervals;  21  January;  (o-p)  and  (q-r) 
1200Z  and  1800Z  respectively,  (s-t)  22/OOOOZ  January  1989. 


20/122  MQI  sip  obs _  20/122  NORAPS  sip  obs 


Figure  20.  (Continued) 
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21/OOZ  MQI  sip  obs  21/OOZ  NORAPS  sip  obs 


21/0 6 Z  MQI  Sip  Obs  2 1/0 6 Z  NO RAPS  sip  obs 


As  shown  in  Fig.21,  the  ridge  line  of  the  mountmn  range  runs  southwest  to  northeast 
along  the  northwest  Maine  border.  From  20/2  lOOZ  to  21/OOOOZ,  the  primary  cyclone 
appeared  to  “jump”  over  the  mountains  to  the  southeastern  side  of  the  ridge  near  44. 8N 
69.5W  (central  Maine)  by  21/OOOOZ  (Fig.20a).  This  indicates  a  speed  of  movement  of 
approximately  25  kts.  The  central  pressure  dropped  2  mb  to  995  mb  by  this  time.  However, 
it  is  not  clear  as  to  whether  this  low  jumped  across  the  mountains  to  combine  with  the  lower 
pressure  area  near  the  coast,  or  dissipated  while  the  secondary  low  pressure  area  developed 
into  the  primary  lOP  5A  cyclone.  The  sparse  data  in  this  mountainous  area  may  make  it 
impossible  to  conclusively  determine  which  scenario  best  fits. 

NORAPS  static  stability  fields  for  the  time  period  20/1 500Z  to  21/OOOOZ  (Fig.  22) 
depict  a  agnificantly  lower  static  stability  over  the  New  England  coastal  region  than  what  was 
over  the  primary  lOP  5A  cyclone’s  path  north  of  New  York,  Vermont  and  New  Hampshire. 
This  reduced  static  stability  along  the  coast  favors  continued  development  of  the  coastal  low 
pressure  area.  By  21/OOOOZ,  NORAPS  moves  the  northern  low  to  the  east  near  46.6N  70.5W 
and  deepens  it  to  996  mb.  The  coastal  low  continued  to  track  along  the  coast  to  near  44.0N 
69.2W.  This  position  places  it  approximately  35  nm  south  of  the  analyzed  low.  A  northwest 
to  southeast  model  cross  section  bisecting  the  secondary  low  along  the  coast  at  21/OOOOZ 
is  shown  in  Fig.  23.  Evident  in  the  cross  section  is  the  strong  low-level  winds  crossing  the 
northern  Appalachian  Mountains  and  the  model  low  pressure  area  located  on  the  leeward  or 
coastal  side  of  the  mountains.  A  66  m/s  jet  streak  (into  the  paper)  is  seen  at  approximately 
300  mb  in  the  top  right  side  of  the  figure.  This  position  places  the  developing  coastal  low 
pressure  tystem  in  a  favorable  position  with  the  transverse  circulation  associated  with  the  jet 
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steak  as  evidenced  by  the  strong  upward  vertical  motion  over  the  low  center  (Uccellini  et. 
al  1985;  Wash  et  al  1988). 

From  these  detailed  comparisons  of  the  analyses  and  forecast  model,  it  is  concluded 
that  lOP  5A  evolved  from  two  separate  low  centers  that  the  NORAPS  model  portrayed  in  an 
accurate  manner.  Lee  side  cyclogenesis  and  strong  upper-level  forcing  enhanced  by  a  lower 
static  stability  likely  played  a  role  in  the  development  of  the  secondary  mesoscale  low  to  the 
south  of  the  primary  cyclone.  This  low  deepened  rapidly  and  likely  became  the  new  lOP  5  A 
cyclone.  The  ability  of  the  high  resolution  NORAPS  model  to  accurately  forecast  the 
secondary  low  confirms  the  speculation  that  previous  failed  attempts  to  model  this  low  was 
due  to  a  lack  of  adequate  resolution. 

C.  EXPLOSIVE  DEVELOPMENT  STAGE  (21/OOOOZ  -  21/1800Z) 

During  this  period,  lOP  5A  experienced  its  greatest  intensification,  deepening  31  mb 
in  18  h.  From  21/OOOZ  to  21/0300Z,  the  primary  cyclone  continued  to  move  to  the  southeast 
to  approximately  44.5®N  67.9°W  and  deepened  4  mb  to  991  mb  (Fig.  20g).  NORAPS 
combined  the  two  previous  low  centers  to  form  the  new  lOP  5A  cyclone  at  the  entrance  to 
the  Bay  of  Fundy  near  44.5  °N,  67.0  °W  (Fig.  20h). 

By  21/0600Z  (Fig.  12g),  the  southern  portion  of  the  Arctic  front  has  dissipated  and 
is  now  located  just  west  of  the  cyclone  center.  There  is  now  a  noticeable  cyclonic  rotation 
of  the  Arctic  front  towards  the  low  center.  The  synoptic  cold  front  continues  to  move  east 
and  the  associated  warm  front  pushes  further  to  the  northeast.  The  previous  warm  front 
intersects  the  Arctic  front  in  New  Brunswick  and  extends  to  the  northeast  into  the  Gulf  of  St. 
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and  warm  fronts.  The  cyclonic  curvature  of  the  Arctic  front  is  also  evident  in  the 
frontogenesis  field  over  northern  Maine  and  New  Brunswick. 

Very  strong  cold-air  advection  accompanied  by  convective  cloud  lines  in  the 
21/0601Z  satellite  imagery  (Fig.  25)  were  present  offshore  from  North  Carolina  to 
Massachusetts  in  the  cold  air  behind  the  front.  This  area  corresponds  well  with  the  large  (> 
1600  watts/m^)  surface  latent  heat  fluxes  simulated  by  NORAPS  in  this  area  (Fig.  26).  Not 
shown  are  the  large  dewpoint  depressions  observed  in  this  area  which  suggests  the  potential 
for  large  latent  heat  fluxes.  The  limiting  streamline  along  the  western  cloud  boundary  has 
continued  to  grow  with  colder  cloud  top  temperatures  evident  over  the  warm  conveyor  belt 
during  the  last  6  h.  A  model  cross  section  (location  shown  in  Fig.  25),  which  bisects  the 
northern  portion  of  the  warm  conveyor  belt  is  shown  in  Fig.  27.  The  stippled  region  are 
model  saturated  values  (>  80%).  NORAPS  does  an  excellent  job  in  depicting  the  dry  edge 
of  the  limiting  streamline.  The  contrast  in  the  moisture  content  of  the  two  air  streams  is 
readily  evident  in  the  gradient  of  relative  humidity  across  the  limiting  streamline,  which 
separates  the  moister  airstream  ahead  of  the  trough  from  the  dryer  air  in  the  upper 
troposphere  west  of  the  trough. 

By  21/1200Z,  the  lOP  5A  cyclone  had  moved  from  the  Bay  of  Fimdy  northeast  to  a 
position  on  the  eastern  side  of  Prince  Edward  Island  in  the  Gulf  of  St.  Lawrence.  A 
deepening  rate  of  13  mb  h'*  from  21/0600Z  to  21/1200Z  is  observed,  and  the  central  pressure 
dropped  to  975  mb  (Fig.l2i).  Two  observations  in  the  area  of  the  MQI  low  report  pressures 
as  low  as  972  mb  (Fig.  20o).  Thus,  the  NORAPS  pressure  of  973  mb  appears  to  be  more 
accurate  (Fig.  12p).  However,  NORAPS  positions  the  low  approximately  60  nm  north  of 
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Figure  23.  Model  cross  section  of  potential  temperature  (solid  contour, 
K),  section  normal  wind  speed  (contoin  interval  20  m  s'‘)  and  cross 
section  circulation  vectors  (m  s"')  along  line  of  Fig.  12f  at  21/OOOOZ 
January  1989. 


Figure  24.  850  mb  temperature  as  in  Fig.  13,  except  for  21/0600Z  January  1989. 


Prince  Edward  Island.  The  Arctic  front  continues  to  intensify  and  move  east  with  a  now 
pronounced  cyclonic  turning  toward  the  low  position.  The  synoptic  cold  front  continues  to 
move  eastward  with  the  occlusion  nearly  intersecting  the  Arctic  front  at  the  low  position. 
The  observations  (Figs.20  m,o)  have  been  showing  colder  temperatures  for  the  last  6  h  ahead 
of  the  occlusion  and  warm  front,  verifying  the  warm  occlusion  classification.  The  surface 
temperature  field  south  and  east  of  the  low  center  is  showing  a  cyclonic  rotation  of  the 
thermal  ridge. 
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Figure  25.  Goes  enhanced  IR  satellite  imagery  at  21/0601Z  January  1989 


The  rotation  of  the  thermal  ridge  is  more  apparent  at  the  850  mb  level  as  depicted  by 
NORAPS  (Fig.  28). 

A  fracture,  or  weakening  of  the  surface  temperature  gradient  along  the  northern 
portion  of  the  cold  front,  similar  to  the  Shapiro  and  Keyser  (1990)  T-bone  model  is  becoming 
evident  in  the  surface  isotherm  pattern.  Similar  weakening  is  also  in  the  surface  frontogenesis 
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Figure  26.  Surface  latent  heat  flux  (contour  interval  100  watts/m^),  valid  at  21/0600Z 
January  1989. 

fields  (Fig.  14e).  Hovv^ever,  the  position  of  the  surface  low  west  of  this  fracturing  indicates 
that  the  cold  front  might  be  catching  up  with  the  warm  front  as  seen  by  Mass  and  Schultz 
(1993),  who  found  similar  weakening  of  the  Lagrangian  frontogenesis  field  in  the  northern 
portion  of  the  cold  front,  but  found  no  similar  fracturing  of  the  surface  thermal  field  for  a 
continental  cyclone  case.  However,  above  the  boundary  layer  (800  mb)  Mass  and  Schultz 
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Figure  27.  Model  vertical  cross  section  of  potential  temperature  (dark 
solid,  contour  interval  2  °K);  relative  humidity  (light  solid,  contour 
interval  20%;  stippled  area  >80%);  and  section  normal  isotachs 
(dashed,  contour  interval  20  ms"')  valid  at  21/0600Z  January  1989. 
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Figure  28.  850  mb  temperature  as  in  Fig.  13,  except  for  21/1200Z  January  1989. 


found  similar  thermal  and  firontogenetical  structures  to  those  of  the  T-  Bone  model  of 
Shapiro  and  Keyser  (1990).  Mass  and  Schultz  (1993)  believe  one  possible  answer  for  the 
differing  surface  structures  is  in  the  profoimd  difference  in  the  trajectories  of  air  parcels 
within  the  cold  and  warm  frontal  zones.  Another  possible  answer  may  be  due  to  the  enhanced 
friction  over  land,  which  contributes  some  frontogenetical  forcing  that  helps  maintain 
continental  cold  fronts  near  the  surface.  However,  over  the  ocean,  especially  over  the  western 
basins  during  winter  where  there  is  a  rapid  modification  of  the  cold  air  masses  from  the  large 
surface  heat  fluxes,  a  frontolytic  environment  is  produced.  Thus,  temperature  gradients 
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within  the  area  of  weak  frontogenesis  along  the  northern  cold  front  can  attenuate  rapidly, 
producing  a  fractured  cold-frontal  surface  (Mass  and  Schultz  1993).  Satellite  imagery  (Fig. 
29)  valid  at  21/1201Z,  shows  a  distinct  mature  comma  head  cloud  with  the  largest  cloud  mass 
still  well  ahead  of  the  low. 

A  vertical  cross  section  at21/1200Z  from  Bedford,  Massachusetts  (BED)  42.4° 
N,  71.3  “W  (Fig.  30)  passes  south  of  the  surface  low,  crosses  the  warm  tongue  near  station 
N1 1 1  then  terminates  at  St.  Johns,  Newfoundland  (YYT)  47.7°  N,  52.8  °W.  The  location  of 
the  cross  section  is  shown  in  Fig.  12j.  Starting  from  the  east,  this  cross  section  verifies  the 
location  of  the  Arctic  front  to  be  between  850  mb  and  800  mb  near  Halifax,  Novia  Scotia 
(YAW).  This  cross  section  continues  to  the  northeast,  intersecting  the  occluded  front  near 
station  Nlll,  approximately  10  nm  southeast  of  Novia  Scotia.  NORAPS  shows  the  Arctic 
front  to  be  somewhat  shallower  between  850  mb  and  925  mb  ^ig.  31).  Otherwise,  both  the 
observed  and  modeled  cross  sections  are  in  excellent  agreement  with  each  other  in  the  thermal 
structure.  NORAPS  moisture  values  show  a  deep  moist  tongue  to  925  mb  east  of  the 
occlusion.  This  structure  is  also  seen  in  the  observed  cross  section.  The  model  wind  fields 
also  correspond  well  with  the  observed  fields  again  noting  the  shift  in  the  winds  over  the 
occlusion. 

The  lOP  5  A  cyclone  continues  to  rapidly  deepen  from  21/1200Z  to  21/1 800Z,  at  a 
rate  of  1 1  mb  6  h'\  This  deepening  rate  is  significantly  greater  than  the  3  mb  6  h‘*  rate  shown 
in  the  NORAPS  coarse  grid  01  analyses  for  this  period.  During  this  time,  the  975  mb  low 
moved  rapidly  northeast  across  the  Gulf  of  St.  Lawrence  at  30  lets  and  deepened  to  964  mb 


Figure  29.  Goes  enhanced  IR  satellite  imagery  at  21/1201Z  January  1989. 


(Figs.  20  o,q).  NORAPS  also  presents  a  westward  bias  in  the  track  with  a  105  nm  error  to 
the  northwest.  There  is  now  a  significant  change  in  the  frontal  evolution.  The  Arctic  front  has 
weakened  and  combined  with  the  previous  baroclinic  zone  along  the  ice  edge  to  the  east. 
The  baroclinic  zone  to  the  east  was  advected  onshore  by  the  strong  southeasterly  flow  caused 
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Figure  30.  Analysis  vertical  cross  section  of  potential  temperature  (dark  solid,  contour 
interval  1  °K);  relative  humidity  (dashed,  contour  interval  10%)  valid  at  21/1200Z  January 
1989. 


by  the  stronger  pressure  gradient  from  the  eastward  moving  cyclone.  There  is  now  strong 
cyclonic  turning  of  the  Arctic  front  into  the  low  center,  giving  the  appearance  of  a  bent-back 
warm  front  in  the  mature  stage.  This  bent-back  thermal  structure  is  evident  at  the  surface 
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Figure  31.  Model  vertical  cross  section  of  potential  temperature  (dark  solid, 
contour  interval  1  °K);  relative  humidity  (light  solid,  contour  interval  20%; 
stippled  area  >80%)  valid  at  21/1200Z  January  1989. 
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but  more  so  at  the  925  mb  (not  shown)  and  850  mb  level  (Fig.32)  Cold  air  encircling  aroimd 
the  cyclone  center  has  now  totally  enclosed  a  region  of  warmer  air  forming  a  warm-core 
seclusion,  similar  to  the  seclusion  process  as  described  by  Shapiro  and  Keyser  (1990).  This 
seclusion  process  is  different  than  Bjerknes  and  Solberg's  (1922)  seclusion  process,  which 
resulted  from  the  topographic  retardation  of  the  front.  Frontogenesis  fields  for  this  time 
depicts  strong  frontogenesis  values  to  the  north  and  southwest  of  the  low  center  correlating 
well  with  the  maximum  baroclinicity  and  wind  speeds  found  in  the  observations  west  of  the 
low  (Fig  20q).  The  previous  stationary  front  along  the  northern  portion  of  the  Gulf  of  St. 
Lawrence  has  dissipated.  Satellite  imagery  at  21/1 701Z  (Fig.  33)  shows  a  mature  cyclone 
vortex  cloud  pattern  with  the  coldest  cloud  top  temperatures  ahead  of  the  low. 

A  vertical  cross  section  from  Bedford,  Massachusetts  (BED)  42.6°N,  71.1  ®W, 
crossing  the  southern  extent  of  the  Arctic  front,  through  the  warm  surface  low  and  front,  and 
terminating  in  Gander,  Newfoundland  (YQX)  49.0 °N,  54.8  °W,  is  presented  in  Fig.  34.  The 
location  of  the  cross  section  is  shown  in  Fig.  12k.  The  model  cross  section  is  shown  in  Fig. 
35.  A  pool  of  Arctic  air  dominates  the  west  side  of  the  cross  section  between  900  mb  and  700 
mb.  There  is  wanning  of  the  theta  air  below  950  mb  as  it  passes  over  the  relatively  warmer 
waters  of  the  Gulf  of  Maine.  The  southern  extent  of  the  Arctic  front  is  evident  between 
Halifax,  Novia  Scotia  (YAW)  and  Queensport,  Novia  Scotia  (WOQ)  where  the  tightest 
packed  theta  lines  intersect  the  surface.  The  surface  warm  tongue  in  Fig.  12  k  is  seen  in  the 
temperature  ridge  between  stations  WOQ,  and  YQX.  There  is  a  diying  out  in  the  eastern  part 
of  the  cross  section  as  evidenced  by  the  sharp  relative  humidity  gradient.  This  structure 
resembles  the  warm  conveyor  belt  structure  described  earlier  in  Fig.27.  The  model  cross 


62 


Figure  32.  850  mb  temperature  as  in  Fig.  13,  except  for  21/1800Z  January  1989. 

section  (Fig.  35)  shows  similar  structure  in  the  vertical.  A  dry  tongue  (<  20%  RH)  extends 
down  to  750  mb  which  separates  the  very  moist  areas  of  the  warm  and  cold  conveyor  belts. 
Below  the  dry  slot  high  relative  humidity  values  (>90%)  dominate  near  the  surface. 

The  Arctic  front  location  is  also  evident  in  the  observed  sounding  taken  at  Halifax, 
Novia  Scotia  (YAW)  as  illustrated  by  the  frontal  inversions  and  deep  moisture  layers  (Fig. 
36).  The  model  sounding  at  this  location  also  depicts  a  low-level  inversion  with  high  moisture 
values  (Fig.  37)  closely  resembling  the  actual  sounding. 
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Figure  33.  Goes  enhanced  IR  satellite  imagery  at  21/1 801Z  January  1989. 


Figure  35.  Model  vertical  cross  section  of  potential  temperature  (dark  solid, 
contour  interval  1  °K);  relative  humidity  (light  solid,  contour  interval  20%; 
stippled  area  >80%)  valid  at  21/1800Z  January  1989. 


D.  MATURE  CYCLONE  STAGE  (21/1800Z  -  22/OOOOZ) 

In  the  final  stages  of  development,  the  lOP  5A  cyclone  moves  east  over 
Newfoundland  to  48.0°N,  55.9°W,  and  fills  one  millibar  to  965  mb  by  22/OOOOZ.  The  Arctic 
fi-ont  and  occluded  fi-ont  are  now  clearly  separated  as  the  cyclone  moves  further  to  the  east. 
The  surface  warm  tongue  has  continued  a  cyclonic  rotation  and  the  system  looks  noticeably 
more  like  a  mature  occluded  cyclone.  However,  inspection  of  the  surface  isotherms  (Fig.  20 
m)  reveals  the  continued  presence  of  the  warm-core  seclusion  noted  at  21/1 800Z.  NORAPS 
shows  this  feature  is  present  at  the  surface  and  850  mb  level  (Fig.  38)  as  well.  NORAPS  has 
also  positioned  the  low  near  St.  Georges  Bay,  130  nm  west  of  the  analyzed  low.  Model 
pressure  is  2  mb  weak  at  967  mb.  The  separation  of  the  Arctic  and  occluded  fi’onts  is  also 
evident  in  the  surface  firontogenesis  pattern  (Fig.  14g).  Frontogenesis  has  noticeably  weaken 
north  of  the  low  where  6  h  earlier  it  was  at  its  strongest  values. 

Satellite  imagery  at  21/2100Z  (Fig.  39)  shows  a  dissipating  comma  head  with  a  front 
still  visible.  The  main  cloud  band  to  the  east  of  the  low  has  moved  away  from  the  cyclone 
center  to  the  east.  The  system,  now  fully  mature  and  vertical  begins  to  fill  as  it  moves  out 
over  the  North  Atlantic  Ocean. 
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Figure  38.  850  mb  temperature  as  in  Fig.  13,  except  for  22/OOOOZ  January  1989. 
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Figure  39.  GOES  enhanced  IR  satellite  imagery  at  21/21 OIZ  January  1989. 


V.  SUMMARY  AND  RECOMMENDATIONS 


A.  SUMMARY 

Using  both  observational  data  and  output  from  the  NORAPS  mesoscale  model,  this 
paper  describes  the  frontal  evolution  and  mesoscale  structure  of  an  intense  cyclone  that 
evolved  from  a  coastal  low  pressure  system  that  eventually  became  the  lOP  5 A  cyclone.  The 
storm  simulation  was  quite  successful,  skillfully  capturing  the  intensity  and  movement  of  the 
cyclone,  as  well  as  most  of  its  structural  characteristics.  Problems  with  the  simulation  include 
modest  position  and  intensity  errors  during  the  latter  portions  of  the  simxilation  as  well  as 
minor  differences  with  the  observed  moisture  structure. 

The  model  simulation  of  the  January  1989  lOP  5  A  cyclone  as  well  as  the  observed 
storm  itself,  suggest  both  similarities  and  differences  with  the  T-bone  conceptual  model  of 
Shapiro  and  Keyser  (1990).  Like  the  Keyser-Shapiro  model,  thermal  seclusion  and  some 
fracturing  of  the  cold  front  were  present  at  the  surface  and  even  more  so  at  850  mb. 
However,  the  storm  also  shows  a  classic  occlusion  extending  to  the  southeast  of  the  storm 
(Fig.  12  i).  In  addition,  the  evolution  of  these  frontal  features  is  quite  different  than  other 
studies.  In  particular,  pre-existing  frontal  zones  intensify  during  the  rapid  cyclogenesis  and 
evolve  into  the  cyclone  frontal  structures  rather  than  forming  as  a  result  of  the  cyclogenesis. 

A  composite  of  the  Arctic  front  evolution  is  presented  in  Fig.  40  at  12  h  intervals 
from  20/1200Z  to  22/OOOOZ.  The  Arctic  front  was  a  well  defined  baroclinic  zone  at  the 
incipient  stages  of  development.  This  front  intensified  as  the  cyclone  developed  and  was 
consistently  located  in  the  immediate  vicinity  of  the  low  throughout  the  lOP.  The  southern 
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portion  of  the  Arctic  front  appeared  to  bend  southeastward  or  wrap  cyclonically  around  the 
low  center  in  the  later  stages  of  development.  The  portion  of  the  polar  front  south  of  the 
incipient  low  weakened  as  it  moved  over  the  warmer  waters  and  joined  with  the  pre-existing 
coastal  front. 

Figure  41  shows  a  composite  of  the  coastal  front  positions  eveiy  6  h  from  20/1200Z 
to  22/OOOOZ.  The  coastal  front  began  as  a  weak  baroclinic  zone  along  the  Mid-Atlantic 
states,  located  well  to  the  southeast  of  the  incipient  low.  The  front  continued  to  rotate 
cyclonically  and  intensify,  in  the  presence  of  a  strong  warm  tongue,  conveyor  belt  and 
developing  low.  The  northern  portion  of  the  front  occluded  and  extended  into  the  low  by 
21/OOOOZ.  The  warm  front  associated  with  the  incipient  system  moved  northeast,  weakened 
and  became  stationary  along  the  strongly  baroclinic  coastal  zone  along  the  southern  border 
of  Quebec. 

Additionally,  ERICA  data  was  used  to  determine  the  mesoscale  structure  of  the  over¬ 
land  coastal  development  period.  The  observed  pressure  and  wind  fields  verified  that  coastal 
cyclogenesis  occurred  prior  to  explosive  cyclogenesis.  The  mesoscale  vortex,  which  appears 
to  have  formed  as  the  result  of  topographic  effects,  developed  south  of  the  primary  cyclone, 
deepened  and  tracked  along  the  Maine  coast  and  likely  became  the  lOP  5 A  cyclone.  Lower 
static  stability,  and  strong  upper-level  forcing  from  a  jet  streak  combined  to  aid  in  the 
development  and  subsequent  intensification  of  the  coastal  system. 
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Figure  41.  Model  composite  synoptic  front  positions  in  six  hour  increments  from  20/1 2Z 
to  22/OOZ  January  1989. 


74 


B.  RECOMMENDATIONS 


The  success  of  the  nested  NORAPS  simulation  encourages  further  numerical 
experimentation  with  other  ERICA  type  storms  to  determine  the  different  types  of  frontal 
evolutions  observed.  Specifically,  further  diagnosis  of  the  development  of  the  warm-core 
seclusion  and  cyclone  center  would  aid  in  determining  the  intensity  and  depth  of  the 
mesoscale  circulations  associated  with  the  warm  core  structure.  Additionally,  trajectories 
amilar  to  those  of  Mass  and  Schultz  (1990)  derived  from  a  model  simulation,  could  explain 
many  of  the  structural  elements  of  the  storm  (i.e.,  conveyor  belts  and  fronts)  by  using  air- 
parcel  trajectories  associated  with  these  features.  Further  diagnostics  of  the  frontogenesis 
function  can  determine  which  components  were  most  important  during  the  break  along  the 
northern  section  of  the  cold  front  and  bent-back  warm  front  structure  in  the  latter  stages  of 
development.  Also  an  adiabatic  mountainless  simulation  could  quanitatively  determine  the 
roles  of  the  topography  or  weaker  stability  along  the  baroclinic  coastal  area  in  the  mesoscale 
cydogenesis.  The  success  of  mesoscale  analyses  and  numerical  forecasts  provides  the  tools 
to  understand  and  predict  future  cases  of  explosive  coastal  development. 
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